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Application of the (3(iV)-hyperspherical harmonics to the study of the 
continuum limits of one-dimensional cr-models and to the generation of 
high-temperature expansions in higher dimensions 
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In this talk we present the exact solution of the most general one-dimensional 0( A)-invariant spin model taking 
values in the sphere S -1 , with nearest-neighbour interactions, and we discuss the possible continuum limits. All 
these results are obtained using a high-temperature expansion in terms of hyperspherical harmonics. Applications 
in higher dimensions of the same technique are then discussed. 



We will report here some new results concern- 
ing the study of one-dimensional u-models [1] 
and the generation of high-temperature series in 
higher dimensions [2]. 

Both studies rely on the use of hyperspheri- 
cal harmonics, for which some new formulae have 
been obtained [1]. Let us first review the defini- 
tion of the hyperspherical harmonics. Define in 



the angular momentum 



L afi = i (x a d,3 - x^da) 
and 

L 2 = J2 L afi L a P . 
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The hyperspherical harmonics are simply the 
eigenfunctions of L 2 on the unit sphere S N_1 C 



L 2 Y lm (a) = A,y /m (<7) 



(3) 



with A, = /(/ + N - 2), a £ 5 Ar - 1 . The index 
/ = 0, 1, . . . labels the eigenspaces Ei, and we 
will refer to it as spin. Each eigenvalue is highly 
degenerate: the index m labels an orthonormal 
basis within each E\. The choice of this basis 
is completely arbitrary but for our purposes we 
have foit und useful to use an overcomplete set in 
terms of Cartesian multipoles, i.e. 



Na c 



(4) 



Y, 



p = VN(N + 2)/2 («7<V - ±6°^ (5) 
= in{a ai ...a°" -Traces) (6) 



Y\ ; a i . . . a i 

where "Traces" is such to make Y/ )ai ... ai 
pletely symmetric and traceless and 

1/2 



2'T(/ + A/2) 



/!T(A/2) 

is a normalization factor. Using this representa- 
tion we have been able to obtain a general expres- 
sion for the Clebsch-Gordan coefficients appear- 
ing in symmetric tensor products of irreducible 
representations and scalars built from them. If 
we define 



a 



dfi(<7)Y /imi (a)Y hm2 (a)Y hms (a) (8) 



(9) 



we have computed exact expressions for 



E 



/"<<i<2<3 fhhh 



m i , m 2 , m 3 



and for the 6-j symbols 1Z{1\, I2, W, h, h, h) when 
one of the spins / 8 - is 1 or 2. 

Using the hyperspherical harmonics one can 
study the possible continuum limits of a generic 
one-dimensional u-model [1]. Consider the most 
general Hamiltonian with nearest-neighbour cou- 
plings: 



H = ^2/V{(T x ■ (T x+ l) 



(10) 
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with a x £ S N 1 . The starting point is the ex- 
pansion of the Boltzmann weight 

(i + f> ; (/3)^K)-^K+i)) (n) 
;=i 



where «;(/?) are coefficients that depend on the 
explicit form of V and satisfy |i>;(/3)| < 1. In 
terms of these quantities it is easy to compute 
the two-point function in the spin-A; channel. We 
get 



(Y k (* ) ■Y k (<T x )) = N k v k (p) 



(12) 



where Af k is the dimension of E k . If < v k (p) < 
1, the two-point function is a pure exponential 
and we define a mass m k (p) as 



m k {p) = - log v k (j3) 



(13) 



To study the critical limit we must investigate 
the limit [3 — ► oo (no critical point can exist for 
finite [3 in one dimension). Thus the problem of 
studying the possible continuum limits is reduced 
to the determination of the asymptotic behaviour 
of v k (p) for [3 — ► oo. If in this limit m k -^C^O 
the corresponding continuum two-point function 
is trivial (white noise): 



(Y k (a ) ■ Y k (a x )) cont = const x 6(x) 



(14) 



Thus a non-trivial continuum limit is obtained 
only if m k (p) — ► (i.e. v k (p) — ► 1) as [3 — ► oo at 
least for some values of k. The various universal- 
ity classes are then characterized by the limiting 
mass ratios m k /mi. 

A detailed study [1] has shown the existence 
of infinitely many families of universality classes. 
The standard N-vector universality class is ob- 
tained, for instance, by considering functions V(t) 
which have a unique absolute minimum at t = 1. 
In this case one can show that 



m k (P) = A(p)X k 



(15) 



with A(/3) —> for [3 — ► oo, so that in the contin- 
uum limit we have 



m k 
mi 



^ k 



(16) 



Analogously it is easy to recover the RP uni- 
versality class: if, for instance, V(t) has absolute 
minima at t = ±1 with V'(±l) ^ 0, then we can 
prove that, in the limit [3 — ► oo, 



(17) 
(18) 
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m 2l 


Ml 


m 2k+ i 




= oo 


m 2l 





Besides these two cases there are many other pos- 
sibilities. If, for instance, V(t) has minima at 
t = to ^ ±1 and t = 1 and some technical condi- 
tions are satisfied [1] (physically they correspond 
to requiring that the contributions of the two min- 
ima to the free energy are comparable in the con- 
tinuum limit), we get new universality classes in 
which the mass ratios are not related to the eigen- 
values \ k of L 2 . 

All these results can be interpreted in another 
framework. In one dimension a continuum field 
theory is simply a continuous-time Markov pro- 
cess on the target manifold. Now, the generator 
of a continuous-time Markov process is the convex 
combination of a diffusion part (a second-order el- 
liptic operator) and a jump part (a positive ker- 
nel) [3]. Physically, this means that the particle 
diffuses for a while according to the specified dif- 
ferential operator, and then, at exponentially dis- 
tributed random times, jumps according to the 
specified probability kernel. On the sphere S N_1 
for N > 3, the only 5'0(#)-invariant second- 
order operator is the Laplace-Beltrami operator 
L 2 (and multiples thereof). The Markov process 
corresponding to pure diffusion on S N_1 is the 
standard TV-vector universality class. One can 
consider also processes with a jump part. It is 
easy to see that there is an infinite- dimensional 
family of possible 5'0(#)-invariant jump kernels 
K: indeed, one can specify an arbitrary probabil- 
ity distribution of jump angles 6 £ [0,7r]. Corre- 
spondingly one finds an infinite- dimensional fam- 
ily of possible continuum limits. For instance the 
case we were mentioning above [V(t) with min- 
ima at t = to and t = 1] corresponds to a process 
with diffusion and a jump distribution 8(6 — Oo), 
cos #o = to- Analogous considerations can be ap- 



plied to RP 



N-l 



The standard RP 



N-l 



sality class corresponds to pure diffusion on the 
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Figure 1. Finite-size-scaling function for the 
spin-2 susceptibility for N = 4. The different 
curves correspond to the TV-vector universality 
class (B = 0), to the RP N_1 universality class 
(B = oo) while 5 = 1 corresponds to one of the 
intermediate universality classes introduced in [4]. 



manifold RP N_1 , while processes with diffusion 
and jumps give rise to new universality classes. 

In two dimensions we argued some time ago 
[4] that there exists a one-parameter family of 
universality classes for 0(7V)-invariant u-models 
taking values in S N_1 , which can be obtained 
by taking appropriate limits in a mixed isovec- 
tor/isotensor model. The same universality 
classes are identified in one dimension: besides 
being obtained as limits of a mixed model, they 
also arise as limits of a Hamiltonian of type (10) 
for suitable V. They correspond to Markov pro- 
cesses on S N_1 with diffusion and jumps of ir. 

For all the Hamiltonians (10), we have more- 
over computed [1] the finite-size-scaling (FSS) 
functions for the susceptibility and the correla- 
tion length in periodic boundary conditions, as 
well as the first correction to them. Interest- 
ingly enough, the qualitative shape of the one- 
dimensional FSS functions is very similar to that 
of their two-dimensional counterpart [4]: see Fig- 
ure 1 for an example. 

Let us now turn to the second subject of this 
talk: the problem of generating high-temperature 



series in dimension d > 2. We will use here the 
general method explained in [5]. Consider a the- 
ory with fields <f> and Hamiltonian 

H = $>(^;/?) + $>(^y;/?) (19) 

x (xy) 

where the second sum is extended over all lat- 
tice links (xy), and assume a high-temperature 
("strong-coupling") expansion of the form 



1 + ^2h n (<f> x ,<i> y ;l3) 



(20) 



with h n (<f> x , 4>y', ft) ~ fi n for /?—*■[). In the present 
implementation of the program one must assume 
an additional property: that the integral 

K e^^h n M x A yi \P)---K k ^ x A yk -,m^) 

vanishes if ^ . n 8 - is odd. This property allows 
us to relate every non- vanishing high-temperature 
graph to one or more lattice random walks (RW). 
However, we must notice that there are inter- 
esting models for which this property does not 
hold: for instance, the RP N_1 and the CP N_1 
(T-models. An additional simplification can be ob- 
tained if one requires an orthogonality property: 



e^^hn^Ay^^hn^Ay^P) = 0(22) 



if m ^ «2- In this case all high-temperature 
graphs can be generated by considering only non- 
reversal RW. We refer to [5] for details on the 
practical implementation. 

For this very general class of theories we have 
computed [2] the two-point functions 

dix) = (h^o,^;?)) (23) 

G 2 (x) = (h 2 (<f>o,<f> x ;P)) (24) 

for various lattices in d = 2,3. The result is writ- 
ten as a sum of abstract graphs which have a pair 
of integers (n,L) associated to each link (in [5] 
these diagrams were called skeleta): the integer n 
indicates the spin associated to the link, while L 
is its length. For instance, to the graph appearing 
in Figure 2 we associate the expression 

Li-l L 2 -l 

[ KMl^l+iA) n K 2 ((t>l,(t>l +l ;p) x 



k = 



h = 



(n 3 ,L 3 ) 
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Figure 2. A simple strong-coupling graph. 




n 



>?+i;/?) 



(25) 



;=o 



with <f>l = <f> 2 = <f> 3 and <f> 1 Li = <f> 2 L2 = <f> 3 Ls . This 
expression must then be integrated over all fields 
with measure d<f> exp[g((/>; /?)]. 

To go further one must specify the fields and 
the specific model and compute the explicit con- 
tribution of each skeleton. For 0(A)-invariant 
(T-models with nearest-neighbour interactions we 
write 



H 



V(a x 



y ' 



(26) 



and expand V(cr x ■ <J y ;f3) as in (11); then we ex- 
press all Green's functions in terms of the pseu- 
docharacter coefficients i>/(/3). What must then 
be computed is the group-theoretic factor of each 
skeleton. We have written a program in MATH- 
EMATICA that does this for an arbitrary skele- 
ton. If the spin associated to each link is small 
(/ < 4), the program is very efficient: the most 
complicated skeleta require 1 min of CPU-time 
on a Sun4/330 with 24 Mbytes of RAM. 

The computation can be simplified if one uses 
the explicit knowledge of C'f ; ; and of the 6- j 
symbols 1Z{1\, I2, W, h, h, h)- Indeed, one can ap- 
ply recursively the rules appearing in Figure 3 
where 



Figure 3. 
triangles. 



Reduction relations for bubbles and 



These two rules allow the complete computation 
of nearly all the graphs which appear in our ex- 
pansions; only for a few of the graphs higher-order 
symbols need to be computed. 

This work was supported in part by INFN and 
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